Solid-Liquid Phase Diagrams for Binary Metallic Alloys: Adjustable Interatomic 

Potentials 
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We develop an approach to determining LJ-EAM potentials for alloys and use these to determine 
the solid-liquid phase diagrams for binary metallic alloys using Kofke's Gibbs-Duhem integration 
technique combined with semigrand canonical Monte Carlo simulations. We demonstrate that it 
is possible to produce a wide-range of experimentally observed binary phase diagrams (with no 
intermetallic phases) by reference to the atomic sizes and cohesive energies of the two elemental 
materials. In some cases, it is useful to employ a single adjustable parameter to adjust the phase 
diagram (we provided a good choice for this free parameter). Next, we perform a systematic in- 
vestigation of the effect of relative atomic sizes and cohesive energies of the elements on the binary 
phase diagrams. We then show that this approach leads to good agreement with several experi- 
mental binary phase diagrams. The main benefit of this approach is not the accurate reproduction 
of experimental phase diagrams, but rather to provide a method by which material properties can 
be continuously changed in simulation studies. This is one of the keys to the use of atomistic 
simulations to understand mechanisms and properties in a manner not available to experiment. 

PACS numbers: 81.30.Bx, 82.60.Lf, 02.70.Uu 
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I. INTRODUCTION 

Atomic-scale simulations have become an indispens- 
able tool for modern materials research. For example, 
molecular dynamics (MD) and Monte Carlo (MC) sim- 
ulations are routinely used to investigate phenomena 
that are not easily accessible via experiment or to in- 
terpret experimental results. 1,2 The fundamental input 
to such simulations is a description of the interactions 
between atoms. While first principles methods accu- 
rately describe atomic bonding through quantum me- 
chanical treatments, they are usually limited to a rel- 
atively small number of atoms. Semi-empirical or em- 
pirical interatomic descriptions are often motivated by 
quantum mechanical ideas but represent different ma- 
terials through parameterization schemes in which the 
constants are fit to experimental (and/or first principles) 
data. While this approach has significant problems when 
the resultant potentials are employed under conditions 
for which they were not fitted, they can provide accurate 
results when applied carefully. Such potentials have the 
advantage that they can be easily used for simulations 
involving a very large number of atoms (currently up to 
10 9 ). Since few interesting materials are pure, we focus 
on potentials for metallic alloys in this paper. In order to 
determine the utility of potentials for alloys, we should 
insure that the potentials lead to the correct phases at 
temperatures and compositions of interest. In this paper, 
we describe the determination of binary phase diagrams 
for a particularly flexible choice of potentials for metallic 
alloys. 

While there has been a long tradition of modeling ma- 
terials using potentials that can be adjusted to represent 



different types of materials, such readily adjustable po- 
tentials are commonly pairwise. As such, they do not 
provide a reasonable description of metals (e.g., surface 
relaxation). Ideally, we will be able to easily adjust a 
potential to give the desired metallic phase diagram. 

Examples of easily adjustable pairwise interatomic po- 
tentials include the Lennard- Jones 3 - (L J) and Morsel po- 
tentials. Because of their simplicity and applicability to 
systems with a wide range of properties, these potentials 
have been widely used, both for elemental and multi- 
component systems^^ However, such potentials are only 
realistic for very simple materials, such as noble gases. In 
other materials, bonding is more complex. For example, 
for metals and alloys, it is well known that many-body 
effects play an important rolei^ For such materials, po- 
tentials of the embedded atom method 9 - (EAM) type are 
widely used. Holian et al. 1 ® proposed an extension of the 
Lennard-Jones potential that allows for many body in- 
teractions, like in EAM type potentials. This idea was 
further pursued by Baske o 11 ' 12 : 13 to treat a broad range of 
metallic systems including alloys resulting in a potential 
known as the Lennard-Jones, Embedded Atom Method 
(LJ-EAM) potential. This potential represents a readily 
adjustable description of atomic interactions in metals 
and metallic alloys. Because of the relatively simple an- 
alytical form of the adjustable LJ-EAM potential, it can 
be used as a description of atomic interactions in systems 
with different, controllable thermodynamics properties. 
As such, it provides a ready means to test the influence 
of thermodynamic properties on materials behavior with- 
out the complexity of first principles approaches or all of 
the oversimplifications inherent in pairwise potentials. 

In this paper, we focus on the systematic determina- 
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tion of binary phase diagrams for LJ-EAM potentials. 
The ability to predict phase diagrams for such adjustable 
interatomic potentials is an important first step in iden- 
tifying a set of atomic interactions required to reproduce 
the requisite phases needed to describe specific phenom- 
ena and systems at non-zero temperature. This is also 
important for determining which thermodynamic prop- 
erties are important in particular type of phenomena 
observed experimentally. We were motivated to pursue 
this study through our own attempts to perform molec- 
ular dynamics simulations of liquid metal embrittlcment 
in alloys*^ There have been several contradictory sug- 
gestions as to what type of thermodynamic behavior is 
necessary for this phenomenon to occur. Therefore, we 
specifically consider the solid- liquid regions of the phase 
diagrams for a wide range of metallic binary alloys. We 
focus on two main parameters in describing the alloys 
- relative atomic size and the strength of the chemical 
bonds. Of course, this is not the first attempt to system- 
atically describe binary phase diagrams from an atom- 
istic view. Earlier attempts have examined the phase dia- 
grams of hard sphere 1 - and Lennard- Jones materialsi 16 i 17 
The methods employed to determine the phase behavior 
range from density functional theory^ to Gibbs-Duhem 
integration methods i 19 i 20 We combine LJ-EAM poten- 
tials, molecular dynamics, and Gibbs-Duhem integration 
methods in a Monte Carlo framework to develop a sys- 
tematic understanding of the relationship between poten- 
tial properties and the solid/liquid phase diagram. 



II. 



POTENTIALS 



In this section, we outline the form of the LJ-EAM 
model used in this work. For a binary material described 
by classical LJ pair potentials, the interatomic potential 
between atoms i and j takes the form: 
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where e SiS and a SiS . are the attractive well depth and the 
diameter for the LJ potential describing the interactions 
between species Sj and Sj (s = A or B). The total energy 
of a binary LJ-EAM system is given by the usual EAM 
form: 
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where F Si (pi) is the embedding energy and <j> Sit 
the pair interaction between atoms i and j separated by a 
distance . The embedded function is commonly chosen 
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where the electron density at the site of atom is 
Pi 
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Here, the dimensionless parameter A s represents the 
strength of the many-body term, the parameter (3 S quan- 
tifies the distance over which the electron density decays 
away from an atom position, and Z\ is the coordination 
number of the reference state (e.g., face centered cubic). 

We can combine the EAM form of the total energy 
appropriate for metals with the convenience of the ad- 
justable LJ pair potential by choosing the pairwise term 
in Eq. ([2]) in such a way that the total energy of the ref- 
erence structure as a function of dilation is described by 
a LJ potential. If we include interactions up to second- 
nearest- neighbors, the like atom pair potential 4>aa{t) 
for species A is given by 
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Here, Z2 is the number of second-nearest-neighbor atoms 
and a is the ratio of the second- to first-nearest-neighbor 
distance. The summation of N terms is carried out until 
(1) 4>AA(r) converges. 

Because this potential was fitted to a LJ-form, it has 
only four adjustable parameters, A a, $a, £a (= €aa), 
and a a (= oaa) for the single component material, A. 
This potential form can be easily extended to multi- 
component systems. For EAM binary alloys, we must 
fix seven functions Pa(t), Pb{t), F A (r), F B {r), (j>AA{r), 
4>bb{t) and 4>ab (t) ■ The first six of these are transferable 
from the two monatomic systems. Following Baskes and 
Stanji^ we can obtain the remaining function, 0ab(?")> 
by fitting to a particular alloy structure. Like them, we 
focus on the ordered LIq compound (after correcting a 
small error in Ref. [l3l ). as described in Appendix A: 

1 r F A (^ ln (r)) + F B (^ l0 (r))" 
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If we set A a = and A B = 0, ^,4,4 and 4>ab take exactly 
the same form as the LJ potential. 

<t>AB(r) is fully determined by Eq. (|5|), except for 
(j>AB{ r )- Therefore, we need two parameters, cab and 
a ab , to describe the cross interaction between species A 
and B (in addition to the 8 parameters needed to de- 
scribe pure A and pure B). The cross-species interaction 
parameters e and a are, therefore, the knobs we use to 
control alloy properties. The Lorentz-Bertholet mixing 
rules^i are widely used for obtaining alloy parameters in 
LJ systems; i.e., o A b = (p A + ob) /2 and e A B = sJzazb- 
However, applying the Lorentz-Bcrthclot mixing rules to 
determine cab and gab results in phase diagrams that 
do not correspond to those observed for metallic alloys 
(i.e., the solid/liquid two phase region is quite square 
with very limited solubility in the solid and the liquid). 
Even though the Lorentz-Berthelot mixing rules may be 
appropriate for LJ potentials, there is no reason to ex- 
pect this to be true for LJ-EAM potentials. This is be- 
cause in the LJ-EAM potentials the pairwise interaction 
represents just part of the bonding and it is the total LJ- 
EAM potential that must reproduce the LJ potential for 
dilation. Therefore, to construct the pairwise interaction 
term in LJ-EAM, we apply the Lorentz-Berthelot mixing 
rule to 4>AB(r) in Eq. (J!}, rather than to <pAB(. r ) itself. 
Doing this, we find the well depth €ab should be 

e-AB = ^ {A A eA + A B e B ^ 

+ yj\{l - A A )e A (l - A B )e B \. (6) 

The derivation of this rule is given in Appendix B. 

Although LJ-EAM potentials are not as widely 
used as the classic LJ potential, the LJ-EAM po- 
tential has already been widely applied in the 
literaturej 1:L ' 12 i 13 ' 22 i 23 ' 24 i 25 The equilibrium structure of 
a LJ crystal at K is face centered cubic (FCC) and the 
melting point depends solely on the cohesive energy and 
the interaction range of the potential^ However, both 
the ground state structure and melting temperature of 
single component LJ-EAM system are strongly depen- 
dent on the two many-body interaction parameters A 
and /3»ia As the strength of many-body interaction A in- 
creases, the melting points of the pure elements decrease 
although the cohesive energy remains unchanged. 11 Ex- 
perimentally, most FCC metals exhibit normalized melt- 
ing temperatures k B T m /Eo in the range from 0.025 to 
0.04, where k B is the Boltzmann constant and Eq the 
cohesive energy of the solid at zero pressure. If we set 
A = 0.7 and $ = 7 in the LJ-EAM potentials for both 
components, pure A and pure B will be FCC at T = 
and the normalized melting points fall within this tem- 
perature range for FCC metals.— 

The melting point of an elemental FCC solid is strictly 
proportional to the well depth e (for any choice of A and 
p). For the reference element A, we set e = 0.6 eV and 
a = 2.5 A with which the melting point T m — 1405 K 
was obtained. However, this choice of e and a is rather 



arbitrary, since only the ratios of well-depths e B / ca and 
diameters <7 B /<ja are relevant in determining the binary 
phase diagrams for LJ-EAM potentials. The potential in- 
teractions were truncated between the second- and third- 
nearest-neighbors such that the relaxed cohesive energy 
is Eq = 6.0e for the FCC reference state [see Eq. (|4a)) ]. 

Table U shows a comparison of the basic properties of 
the LJ-EAM material with those of several FCC noble 
metals. For the LJ-EAM parameters employed in Ta- 
ble HI the LJ-EAM potential yields reasonable properties 
for the FCC metals except for the bulk modulus. Unlike 
other properties, the bulk modulus is solely determined 
by the second derivative of the energy as a function of lat- 
tice dilation (regardless of the many-body interactions). 
Since the LJ-EAM potential was fitted to a LJ-form (by 
definition), both the LJ and LJ-EAM models yield sim- 
ilar bulk moduli. The bulk modulus of LJ systems is 
known to be too large compared with the FCC metals.— 
One could obtain better bulk modulus values by fitting 
the pair interactions to another form, e.g., the universal 
binding energy relation? 2 ^ 

III. GIBBS-DUHEM INTEGRATION METHOD 

There are several computational methods that can be 
used to determine the phase diagram of a system de- 
scribed by any choice for the interatomic potentials ! 19 > 28 
Here, we directly construct the phase diagram using 
the Gibbs-Duhcm integration technique proposed by 
Kofke . 19 i 20 In this method, two or more coexisting phases 
are simulated (semi-grand canonical Monte Carlo simula- 
tions) independently at the same temperature and pres- 
sure. Once a single point of the coexistence curve be- 
tween two phases is known, the rest of the curve can be 
computed (without any free energy calculations) by in- 
tegrating the equivalent of the Clausius-Clapeyron equa- 
tion for coexistence during the course of the simulations. 
The Clapeyron equation for equilibrium between two bi- 
nary phases (e.g., liquid and solid) at constant pressure 
is given by 

d(3 _ ( X B ~ X b) fn\ 

dt B ~ £ B (i-z B )(w-H»y {<> 

where (3 is the reciprocal temperature, l/k B T, and T the 
absolute temperature, £ B the fugacity fraction of species 
B, £b = Jb/Ja, fi the fugacity of species i in solution, 
x B the mole fraction of species B, and H is the molar en- 
thalpy. The right-hand side of Eq. ([7|) can be integrated 
numerically to find an equation for j3 as a function of £ B 
if we have an initial condition describing the tempera- 
ture, fugacity fraction, enthalpies, and compositions at 
one coexistence point. 

The Gibbs-Duhcm integration method is a more effi- 
cient approach to determining phase equilibrium in solid 
systems than the Gibbs-ensemble method* 2 ^ because 
the important Monte Carlo move is changing the elemen- 
tal identity of a particle rather than inserting or removing 
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TABLE I: A comparison of several key properties for the LJ-EAM potential model (for A — 0.7, (3 — 7) and several FCC 
metals. All properties are normalized by using E (= 6e), r (= s/2a), and Q(= r%/V2). 



property 


normalized quantity 


LJ 


LJ-EAM 


Cu 


Ag 


Au 


FCC metals 


bulk modulus 


BQ/Eo 


8.5 


8.5 


3.0 


3.92 


5.05 


2.0-6.0 


Melting point 


kBT m /E 


0.082 


0.033 


0.033 


0.037 


0.029 


0.024-0.041 


Cauchy discrepancy 


(C12 - C4A)/B 





0.48 


0.30 


0.425 


0.69 


-0.4-0.7 


vacancy formation energy 


e£/e 


1.0 


0.31 


0.37 


0.39 


0.23 


0.2-0.5 


(100) surface energy 


E(ioo)ro/Eo 


0.67 


0.26 


0.24 


0.23 


0.22 


0.15-0.25 



a particle from the system (inserting a particle is a low ac- 
ceptance probability event). The Gibbs-Duhem integra- 
tion approach has been used to tackle a number of mul- 
ticomponent, multiphase equilibrium problems including 
calculating the phase diagram of a binary Lennard-Jones 
fluids and calculating phase diagrams for colloids in 
polymer solutions*^ 

In this paper, we determine the initial coexistence 
point by performing a microcanonical ensemble molecu- 
lar dynamics simulation of an elemental system contain- 
ing both a solid and liquid. 29 Using this approach, we 
directly measure the melting temperature of pure A or B 
(the solid and liquid fractions in the simulation cell evolve 
until the system reaches the equilibrium melting point). 
Using this data as the starting point, we employ Kofke's 
Gibbs-Duhem integration technique 2 -^ within the frame- 
work of semigrand canonical Monte Carlo simulations. 
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IV. GENERIC BEHAVIOR OF BINARY PHASE 
DIAGRAMS FOR LJ-EAM MATERIALS 

In this section, we investigate solid/liquid phase 
diagrams for binary LJ-EAM systems. All of the 
temperature-composition phase diagrams reported below 
are for binary alloys at atmospheric pressure. In particu- 
lar, we focus on (a) cases where the melting points of the 
two elemental systems are identical and (b) cases in which 
the elemental systems have very different melting points. 
The former (similar melting points) is a relatively com- 
mon case while the latter is more rare (although common 
for systems exhibiting liquid metal embrittlement J^) In 
each case, we examine how variations of the ratio of 
Lennard-Jones diameters obI^a and well-depth cb/^a 
affect the phase diagrams. We also compare major trends 
observed in the simulated phase diagrams with those 
measured experimentally (in order to evaluate the appro- 
priateness of the LJ-EAM alloy model to mimic behavior 
in real metallic systems). 

A. Equal melting points 

Figure Q] shows three types of binary phase diagrams 
commonly observed in metallic alloys; we refer to these 
as azeotropes [Fig.QJa)], simple eutectics [Fig.[IJb)], and 
liquid phase miscibility gap systems [Fig. []Jc)] of which 



FIG. 1: Schematic phase diagrams for metallic binary alloys 
of comparable melting points: (a) azeotrope, (b) simple eu- 
tectic, and (c) a combination of a eutectic, a monotectic and 
a liquid phase miscibility gap. The symbols are as follows: L 
refers to a liquid solution of A and B, S to a solid solution of 
A and B, Sa to a solid solution rich in A, and Sb to a solid 
solution rich in B. 



there are several types. We explicitly omit consideration 
of systems exhibiting compounds (for now). It is the 
atomic size mismatch and the cross-species pair interac- 
tion that determine which phase diagram type pertains. 
To investigate the effect of both atomic size mismatch 
and cross-species pair interactions, we considered three 
series of cases: (1) like bonding (e^s = £a = (b) with 
different atomic size ratios ctb/oa, (2) equal atomic sizes 
(erg = a a = <Jab) and equal well-depths for the elements 
(f-B = £a) but vary the well-depth describing the AB in- 
teractions (-ab > and (3) vary both ob/oa and cabI^a- 
The other parameters, A and j3, are fixed at reference 
values 0.7 and 7 as mentioned in Sec. [TT] 

Figure [5] shows temperature-composition phase dia- 
grams for a series of binary systems of type (1) (variation 
of atomic size) . When the atomic size difference is small 
(less than 8%), the solid region of the phase diagram is 
a solution (for the temperature range examined). How- 
ever, as the atomic size difference increases, the degree 
of phase separation increases and the solid forms a misci- 
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FIG. 2: Solid-liquid phase diagrams with various atomic size 
differences. Model alloys have the same melting point with 
zero heat of mixing (es/e^ = 1.0, e^b/ea = 1.0). Solid-liquid 
phase diagrams were calculated with variation of atomic size 
difference (ovb/cta = 1.05, 1.08, and 1.1). In this and sub- 
sequent phase diagrams, we estimate that each data point in 
the phase diagram has error bars of w 10 ° in temperature and 
~ 3% (although these errors have some variation depending 
on the relative stability of the different phases at each point 
along the phase boundaries). 



bility gap that leads to a eutectic phase diagram starting 
at an atomic size difference between 8% and 10%. This 
cross-over from solid solution to eutectic phase diagrams 
occurs at a size difference that is much small than re- 
ported for LJii (14 ~ 15%) or hard sphere systems^ 
(12.5%). Moreover, the range of size difference for this 
type of simple eutectic phase diagram is very narrow; fur- 
ther increase in this difference leads to a miscibility gap 
in the liquid. When a miscibility gap exists in the liq- 
uid, the phase diagram becomes more complex (similar 
to Fig. [He)), with negligible solubility in solid phases. 
Therefore, eutectic phase diagrams with deep eutectic 
point cannot be obtained by increasing size difference in 
LJ-EAM model. This is not the case for LJ or hard 
sphere materials. 

Figure [3] shows temperature-composition phase di- 
agrams for binary mixtures with same atomic size 
(<tb/(7a = 1) and different values of cross-species inter- 
action parameter. When cab/^a > 1, heat of mixing is 
negative (exothermic solid solution) and A and B atoms 
"like each other" . In these cases, the phase diagrams 
form continuous solid solutions over the whole composi- 
tion range and the liquidus appears to be parabolic. The 
maximum temperature for which the solid and liquid co- 
exist increases as cab/^a increases (the heat of mixing 
is more negative the larger the cab/^a ratio). For a bi- 
nary mixture with well-depth ratio of unity cab/^a = 1, 
all of the atoms are indistinguishable - hence, the phase 
diagram would simply be a horizontal line at T = 1405 
K. On the other hand, when cab/^a < 1, the heat of 
mixing is positive (endothermic solid solution) and A- 
B atom "dislike each other". Under this circumstance, 




0.4 0.6 o.s 
Fraction of B 

FIG. 3: Solid-liquid phase diagrams for different values of 
cab. The pure metals, A and B, have the same melting point 
and the same atomic size (sb/^a = 1.0, ub/^a = 1.0). The 
solid-liquid phase diagrams were calculated with different val- 
ues of the cross-species interaction parameter (€ab/ca = 1-1, 
1.05, and 0.95). 



the liquidus curve is concave for 1 < cab/^A < 0.95 (not 
shown). For €ab < 0.95, the two species are no longer 
miscible at all compositions and a large miscibility gap 
appears in both the solid and liquid phase regions, as 
shown in the phase diagram in Fig. [3) 

These results suggest that the heat of mixing alone 
does not control the melting point (liquidus). Rather, 
atomic size difference also plays an important role, es- 
pecially for forming eutectic phase diagrams. However, 
both the size difference and heat of mixing affect the ten- 
dency towards mixing and can act to compensate each 
other. Therefore, controlling atomic size difference to- 
gether with the heat of mixing can lead to a wide range 
of types of binary phase diagrams. For example, increas- 
ing the liquidus temperature by increasing cab/^a can 
be compensated by increasing the atomic size difference. 
When large atomic size difference are combined with rel- 
atively large cab / f-A , eutectic phase diagrams with deep 
eutectics can be formed, as shown in Fig. 2J This trend 
is quite interesting because it seems to be related with 
the rule of thumb for making metallic glasses. In metallic 
glass systems, it has been established, empirically, that 
the ability to form glasses is greatest in multicomponent 
systems in which the atomic size difference is large and 
the heat of mixing is strongly negative (large cab / &a) *~ 
The present phase diagram results suggest that this is 
also the description of the condition for the formation of 
deep eutectics. Deep eutectics are also known as systems 
for which glass formation is particularly easy. 



B. Phase diagrams of large melting point difference 

Solid-liquid coexistence is a key to many materials pro- 
cessing strategies involving solidification and in-service 
conditions where a solid metal is in contact with an- 
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FIG. 4: Solid-liquid phase diagram for a model binary alloys 
with atomic size differences of 10% - 13% and e^s /(J a = 1.03. 
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FIG. 5: Schematic phase diagrams of binary alloys where the 
two components have very different melting points. Two such 
cases are (a) a eutectic with significant solubility in the solid 
and (b) a diagram with two eutectic points and a miscibility 
gap in the liquid phase. 



other, liquid metal. In many of the latter cases, the liquid 
phases consist of low melting point species, such as Hg, 
Ga, Bi, Pb, and Sn. Metallic binary systems in which the 
melting points of the two components differ greatly typ- 
ically show one of two types of simple solid/liquid phase 
diagram (provided no intermetallic compounds form): 
these are eutectics phase diagrams with or without a liq- 
uid phase miscibility gap, as shown in Fig. [5] When there 
is a liquid phase miscibility gap, the solubility of the mi- 
nority species in the solid phase is usually very small [al- 
though it often appears exaggerated in schematic phase 
diagrams, such as Fig. Ob)]. However, in other binary 
systems, such as Al-Ga and Zn-Ga, the solubility in the 
solid phase can be significant over the entire tempera- 
ture range, in spite of large melting point difference [see 
Fig. Ha)]. 

Phase diagrams were calculated for several different 
binary systems, where we varied the melting point of 
species B. Figure [5] shows the temperature-composition 
phase diagrams for systems with the atom size ratio fixed 
as ob/va =1 and well-depth ratios of es/^A — 0.6, 0.5, 
and 0.4. When the melting point of B is comparable to A, 
the system forms a solid solution with a spindle shaped 
solid- liquid two-phase region. As es/^A decreases, A-B 
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FIG. 6: Solid-liquid phase diagrams for binary alloys with 
different melting points of component B, as controlled by dif- 
ferent choices of es/eA(=0.6, 0.5, and 0.4). These alloys have 
the same atomic sizes {<jb/o'a = 1.0, and the parameter sab 
is described by Eq. ((U). 
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FIG. 7: Solid-liquid phase diagrams with different atomic 
sizes ratios (o\b/o"a =0.9, 1.0, and 1.1). The melting point of 
B was held constant by fixing £b/ca = 0.5 and Ab = 0.7). 



become weaker and the phase diagram evolves to a eutec- 
tic diagram. It is difficult to see the solid-liquid two-phase 
region at the B-rich side of the phase diagram in Fig. [5] 
because the eutectic point is close to pure B. Nonethe- 
less, we assure the reader that this is indeed a eutectic, 
just like in Fig.OJa). As €b/^a decreases further, the sol- 
ubility of B in the solid phase and the solubility of A in 
the liquid phase decrease. This trend agrees with the ob- 
servation that eutectic phase diagrams determined from 
experiment tend to show smaller solubilities as the ratio 
of the melting points of the two species deviates further 
from unity. 

Phase diagrams were also determined for binary sys- 
tems with different atomic size ratios. Figure [7] shows 
the temperature-composition phase diagrams for binary 
mixtures with a fixed well-depth ratio, cb/ca — 0.5, 
and several diameter ratios, ub/o'a = 0.9, 1.0, and 1.1. 
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0.4 0.6 0.8 

Fraction of B 

FIG. 8: Solid-liquid phase diagrams with different atomic 
sizes ratios (<tb/<7/i =0.9, f.0, and f.f). The melting point of 
B was held constant by fixing Ab = 0.9 and es/e^ = I. 



When there is no atomic size mismatch, the solubilities 
in the solid and liquid phases are quite large, despite the 
large melting point difference. But, as the atomic size 
mismatch decreases, these solubilities decrease. Interest- 
ingly, when the atomic size of B is larger than that of A, 
the size effect is dominant (i.e., the solubility is negligible 
over the entire temperature range and a miscibility gap 
appears in liquid phase). 

The melting point of B can be set by the choice of 
the many-body interaction parameter Ab (in addition to 
choosing the well-depth ratio €b/za)- The melting point 
of B decreases with increasing A b even though the cohe- 
sive energy remains unchanged. Figure [5] shows the solid- 
liquid phase diagram with Ab = 0.9 for several different 
atomic size ratios. Effect of atomic size mismatch is still 
valid for this kind of solid/liquid pairs. The trends in the 
phase diagrams with atomic size mismatch are similar in 
this case to those shown in Fig. [7l 



TABLE II: Cohesive energy Eo, lattice parameter ao, and 
melting point T m for column IB FCC metal a 31 ' 32 and the cor- 
responding potential parameters e, a, and A. 



element 


Eo (eV) 


ao (A) 


T m (K) 


c 


(7 


A 


Cu 


3.54 


3.62 


1357 


0.59 


2.277 


0.7 


Ag 


2.85 


4.09 


1235 


0.475 


2.573 


0.66 


Au 


3.93 


4.08 


1392 


0.655 


2.570 


0.8 
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C. Comparison with real binary alloy systems 

Since the phase diagrams shown above were deter- 
mined within the framework of generic interatomic po- 
tentials LJ-EAM, it is interesting to inquire to what de- 
gree choosing parameters in the potential can lead to 
phase diagrams that are consistent with those found ex- 
perimentally in real metallic systems. We can compare 
the simulation data to experimental results to verify the 
ability of the LJ-EAM model to mimic behavior in real 
metallic systems. The elements found in column IB of the 
periodic table, copper, silver, and gold, are common FCC 
metals that are well described with the embedded atom 
method framework. We determine the LJ-EAM parame- 
ters, e, a and A, to reproduce the cohesive energy, lattice 
parameter and melting temperatures of these column IB 
elements (see Ref . Il2h , as shown in Table [TTJ 

We employ these parameters within the LJ-EAM 



1100 - 




FIG. 9: Solid-liquid phase diagrams for (a) Ag-Cu, (b) Cu- 
Au, and (c) Ag-Au. Cohesive energy and LJ diameter were 
fitted to real properties of materials and cab/cab was ad- 
justed in LJ-EAM binary alloy model. £ab/(°ab = 1-03 for 
Ag-Cu, e AB /e AB = 1.03' for Cu-Au, e AB /e° AB = 1.0 for Ag- 
Au, where e% is given in Eq. ((6]). 
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framework to calculate the corresponding Ag-Cu, Cu-Au, 
and Ag-Au binary phase diagrams. By only adjusting 
the ratio cab/za [note, the final values were very close 
to those predicted by Eq. (JBJ], we are able to reasonably 
reproduce the experimental phase diagrams (each with 
a unique topology). Although the agreement certainly 
is not perfect, the phase diagram type, the temperature 
ranges of the features of the diagrams, and solubilities 
are in good agreement with experiment. This is remark- 
able given that only one parameter was varied (and the 
atomic sizes and cohesive energies are available from ex- 
periment). This type of agreement is possible for many 
binary systems, provided that they do not exhibit inter- 
metallic compounds. 



V. CONCLUSION 

We developed an approach to determining LJ-EAM 
potentials for alloys and used these to determine the 
solid-liquid phase diagrams for binary metallic alloys 
using Kofke's Gibbs-Duhem integration technique com- 
bined with semigrand canonical Monte Carlo simulations. 
Inclusion of many-body interactions led to phase dia- 
grams which can be quite different from those determined 
using LJ or hard sphere materials. We demonstrated that 
it is possible to produce a wide-range of experimentally 
observed binary phase diagrams (with no intermetallic 
phases) by reference to the atomic sizes and cohesive en- 
ergies of the two elemental materials and by judicious 
choice of a single parameter that controls the pairwise 
interactions of these two elements. In addition, we pro- 
vided a formula that leads to good choices for this one 
free parameter. Within this framework, we performed a 
systematic investigation of the effect of relative atomic 
sizes and cohesive energies of the elements on the bi- 
nary phase diagrams. Finally, we demonstrated that this 
approach leads to good agreement with several experi- 
mental binary phase diagrams. The main benefit of this 
approach is not, in our opinion, to accurately reproduce 
the phase diagrams of real materials. Rather, it is to 
provide a method by which material properties can be 
continuously changed in simulation studies in order to 
develop understanding of mechanisms and properties in 
a manner not available to experiment. To this end, the 
relationship between the adjustable potentials and the 
phase diagrams they imply is central. 
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APPENDIX A: DERIVATION OF <j} AB (r) 

To derive formulation for 4>(r), Llo structure with 
c/a=l was considered as reference state. In LJ-EAM 
formalism, the energy per atom of this structure as a 
function of dilation is given by 



E h L t EAM (r) 



Llo 



+ [(f>AA(r) + Hab{t) + (f>BB{r)} 



+ 2 [4>AA(ar) + (f> B B(ar)] , 



where 



p L A l0 (r) = i [2p B (r) + PA (r)} + \ PA {ar) 



and 



P B °(r) = I i 2 PA(r) + Pb(t)} + l - PB {ar). 



(Al) 



(A2) 



(A3) 



If we rewrite the energy of this system in terms of LJ 
pair potentials, 

E L L i(r)= ^ A (r)+4cb L / B (r)+^ B J B (r). (A4) 

By setting ££f o EAM (r) = E\l a {r), we can derive Eq. 

APPENDIX B: MIXING RULE FOR e AB 

In this appendix, we show the derivation of the ex- 
pression we employ for describing the interaction between 
unlike species. Given Eq. this reduces to the determi- 
nation of 4>V B ■ The length scale parameter a AB is simply 
the arithmetic average of those for the elements. There- 
fore, we only require a mixing rule for the well depth 
parameter €ab- In this appendix, we show the origin of 
the choice for this parameter that was quoted in Eq. ([6]). 

In order to keep this analysis simple, we rewrite 
Eq. (|A1[) for the special case of nearest-neighbor inter- 
actions, a a = ob and (3 a — Pb- 

Eli (t)= \[FA{p\{r))+FB{PB{r))} 

+ [4>AA{r) + 4<p A B(r) + <Pbb(t)} 
= 0&(r)+4^(r)+#&(r). (Bl) 

In equilibrium at zero pressure, r = r eq (= \/2<ta = 
V2<tb), p° A (r eq ) = p° B (r eq ) = 1, ^ A (r eq ) = e A , 
< f > BB( r eq) = e B, and 4> AB (r eq ) = €ab ■ Substituting these 
expressions into Eqs. (|3"a|) and (|4"a[) . yields 



F A {f A {re q )) 



F B {p%{r eq )) 



A A Zie A = 6A A eA, 
A B Z\e B = QA B e B , (B2) 
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and 



4>AA{req) = 0AA.( r eq) , 

, Z 1 



i^(p>e 9 )) 
= £A - A4<U = ( 1 - A.Aj e A , 
<pB B (re q ) =(l-A B )e B . (B3) 
The energy per atom at r=r eq then becomes 



E Li (r eq ) 



- 6A A e A + 6A B e B 

+ [(1 - A A )e A + A(t> AB (r) + (1 - A B )e B 
= e A + 4e AB + e B . (B4) 

Solving this for e AB yields 

e-AB = ^ (A A e A + A B e B ^j + <f> AB {r eq ) (B5) 



Now, determining e AB is a matter of determining 
4>AB{r e q)- As described in the text, we obtain <pAB(r eq ) 
by applying the Lorentz-Bertholet mixing rule to </>(r eg ) 
rather than to e (these are equivalent for the pairwise 
potentials): 

4>AB{r eq ) = \J\<pAA(r eq )^BB(r e q)\ 



(1 - A A )e A (l - A B )e B \. (B6) 



Substituting this into Eq. (|B5|) gives the modified mixing 
rule, Eq. ©: 



A A e A + A B e B 



{l-A A )e A (l-A B )e B \. (B7) 



* Electronic address: hnam@princeton.edu 

1 K. Ohno, K. Esfarjani, and Y. Kawazoe, Computational 
Materials Science (Springer, Berlin, 1999). 

2 D. Frenkel and B. Smit, Understanding Molecular Simula- 
tion (Academic Press, San Diego, 2002), 2nd ed. 

3 J. E. Lennard- Jones, Proc. R. Soc. Lond. A 106, 463 
(1924). 

4 P. M. Morse, Phys. Rev. 34, 57 (1929). 

5 E. J. Jensen, W. D. Kristensen, and R. M. J. Cotterill, 
Philos. Mag. 27, 623 (1973). 

6 J. Q. Broughton and G. H. Gilmer, J. Chem. Phys 79, 
5095 (1983). 

7 V. Vitek and D. J. Srolovitz, eds., Atomistic Simulation 
of Materials Beyond Pair Potential (Plenum, New York, 
1989). 

8 A. E. Carlsson, Beyond Pair Potentials in Elemental 
Transition Metals and Semiconductors (Academic Press, 
Boston, 1990). 

9 M. S. Daw and M. I. Baskes, Phys. Rev. Lett. 50, 1285 
(1983). 

10 B. L. Holian, A. F. Voter, N. J. Wagner, R. J. Ravelo, S. P. 
Chen, W. G. Hoover, C. G. Hoover, J. E. Hammerberg, 
and T. D. Dontje, Phys. Rev. A 43, 2655 (1991). 

11 M. I. Baskes, Phys. Rev. Lett. 83, 2592 (1999). 

12 S. G. Srinivasan and M. I. Baskes, Proc. R. Soc. Lond. A 
460, 1649 (2004). 

13 M. I. Baskes and M. Stan, Metall. Mater. Trans. A 34, 435 
(2003). 

14 B. Joseph, M. Picat, and F. Barbier, Eur. Phys. J. Appl. 
Phys. 5, 19 (1999). 

15 W. G. T. Kranendonk and D. Frenkel, Mol. Phys. 72, 679 



(1991). 

16 M. J. Vlot, J. C. van Miltenburg, H. A. J. Oonk, and J. P. 
van der Eerden, J. Chem. Phys 107, 10102 (1997). 

17 M. R. Hitchcock and C. K. Hall, J. Chem. Phys 110, 11433 
(1999). 

18 W. G. T. Octoby, Mol. Phys. 72, 679 (1991). 

19 D. A. Kofke, Mol. Phys. 78, 1331 (1993). 

20 D. A. Kofke, J. Chem. Phys 98, 4149 (1993). 

21 J. S. Rowlinson, Liquids and Liquid Mixtures (Butterworth 
Scientific, London, 1982). 

22 E. Yurtsever and F.Calvo, Phys. Rev. B 62, 9977 (2000). 

23 T. Wang, F. X. Zhou, and Y. W. Liu, Chinese Physics 11, 
139 (2002). 

24 M. I. Baskes, JOM 56, 45 (2004). 

25 C.-W. Pao and D. J. Srolovitz, Phys. Rev. Lett. 96, 186103 
(2006). 

26 J. A. Pryde, The Liquid State (Hutchinson University Li- 
brary, London, 1969). 

27 J. H. Rose, J. R. Smith, F. Guinea, and J. Ferrante, Phys. 
Rev. B 29, 2963 (1984). 

28 A. Z. Panggiotopoulos, N. Quirke, M. Stapleton, and D. J. 
Tildesley, Mol. Phys. 63, 527 (1988). 

29 J. R. Morris, C. Z. Wang, K. M. Ho, and C. T. Chan, Phys. 
Rev. B 49, 3109 (1994). 

30 A. Inoue, Acta Mater. 48, 279 (2000). 

31 C. Kittel, Introduction to Solid State Physics (Wiley, New 
York, 1995), 7th ed. 

32 C. J. Smith, ed., Metal Reference Book (Butterworths, 
London, 1976), 5th ed. 



